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We apply a type of background independent "polymer" quantization to a free scalar field in a 
flat spacetime. Using semi-classical states, we find an effective wave equation that is both nonlinear 
and Lorentz invariance violating. We solve this equation perturbatively for several cases of physical 
interest, and show that polymer corrections to solutions of the Klein-Gordon equation depend on 
the amplitude of the field. This leads to an effective dispersion relation that depends on the ampli- 
tude, frequency and shape of the wave-packet, and is hence distinct from other modified dispersion 
relations found in the literature. We also demonstrate that polymer effects tend to accumulate 
with time for plane-symmetric waveforms. We conclude by discussing the possibility of measuring 
' deviations from the Klein-Gordon equation in particle accelerators or astrophysical observations. 
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I. INTRODUCTION 



There has been much interest in the past few years on possible quantum gravitational effects on wave propagation. 
The basic idea is that if the metric is quantized, effective wave dynamics follow from replacing the classical metric 
in the wave equation by the expectation value of a metric operator in a suitable quantum state. Such a state would 
provide a "quantum gravity corrected" background spacetime with a built in fundamental discreteness scale. It has 
been conjectured that this may lead to Lorentz symmetry violating modifications to conventional wave propagation 
and dispersion relations Phenomenological effects have been studied by introducing explicit Lorentz violating 
terms in particle physics models with potentially significant effects 0] that could require fine tuning to agree with 
some observations. 

Several authors have previously considered the effects of dispersion relations incorporating high energy corrections 
on Hawking radiation [j, [5[ and inflationary cosmology @. In the studies, the precise form of the corrections was 
specified on an ad hoc basis. It was found that the Hawking effect is not affected by such modifications, whereas in 
f^S . cosmology some classes of dispersion relations do affect the spectrum of quantum fluctuations. 

It has been argued that quantum gravity modifications to the group velocity of wave-packets can lead to potentially 
observable effects on the signal from gamma ray bursts Q ■ Observed violations of the GZK cutoff @ in the spectrum 
of high energy cosmic rays can also indicate a break-down of Lorentz symmetry at high energy. Attempts have been 
Q\ [ made to systematically derive constraints on parameters appearing in modified dispersion relations such as 



E 2 = m 2 + p 2 



1 + C 



i-Y 



(i) 



where E, m and p denote the energy, mass and momentum respectively; p p \ is the Planck momentum and C and 7 > 



are real parameters Q. 

It is important to ask which ideas associated with the the construction of a quantum theory of gravity lead to 
potentially observable effects on wave propagation in the semi-classical regime. The work in this direction utilizes a 
background independent quantization scheme for the dynamical variables associated with geometry. This approach 
to quantum theory may be applied to any system, and in particular also to matter degrees of freedom. Indeed for a 
complete quantization of a gravity-matter system it is natural to apply the same procedure to all degrees of freedom. 

Our purpose in this paper is to study the effects of a background independent quantization on wave propagation 
without any quantization of geometry. The approach may be viewed as the study of polymer quantum fields on curved 
spacetime, which is a general area of interest in its own right. 

The simplest case is the massless scalar field on a Minkowski background. The kinematics of this quantization, 
and its comparison to the usual Fock procedure, has been discussed at a mathematical level [IcJ HH, EH using a 
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direct application of the ideas developed in loop quantum gravity (LQG). However, so far there appears to have been 
relatively little study of possible physical consequences of this approach. This is what we undertake here, but with a 
number of essential differences in approach on which we will elaborate. 

This paper is organized as follows. In fllTl we review some aspects of background independent quantization, with 
emphasis on a new set of basic variables for quantization. (These have application well beyond the scenario addressed 
here.) We then introduce a class of semi-classical states, and show how these lead to a modified wave equation that is 
both nonlinear and Lorentz invariance violating. This has obvious and significant consequences for wave propagation, 
even without considering any quantum gravity effects. It provides a sharp contrast to the linear modifications of the 
wave equation used in all the previous applications mentioned above. We then describe a perturbative solution scheme 
for the semi-classical wave equation in jjllll and its specialization to situations with planar and spherical symmetry. 
In QYV\ we discuss the phenomenological implications of the model in the context of modified dispersion relations, the 
coherence length of collimated beams (with application to the Large Hadron Collider and high energy cosmic rays), 
and pseudo-spherical radiation from astrophysical sources (such as gamma-ray bursts). A summary and discussion of 
our results is in [jV] 



II. QUANTIZATION AND THE EFFECTIVE WAVE EQUATION 



A. Polymer formalism 



We consider a massless free scalar field in four spacetime dimensions with the following background metric 

g^dx^dx" = -dt 2 + q ab dx a dx b . (2) 

Here g^ is the spacetime metric whereas q a b denotes the spatial metric. The canonical phase space variables are 
(</>, P^) with scalar matter Hamiltonian 

H = Jd 3 x (^-LP| + ^Lq^daW^ , (3) 

where q = det{q a b). The usual wave equation for scalar field follows from this via Hamilton's equations. 
Let us consider the functions 

<t>f = J d 3 yVqf(yMy), UxiP^) = exp [iXP^/^j } , (4) 

as the new set of basic variables that are to be realized as operators in the quantum theory. Here f(x) is a scalar and 
A is a real constant with the dimension of length squared in natural units. The factors of ^/q are necessary to balance 
density weights in the integral and in the exponent. (These variables may be viewed as trie "dual" [13fl to those used 
in the polymer quantization of a particle system [ll| [lj| motivated by loop quantum gravity LQG [l4l Il5l . liH llT\ , 
where the exponentiated configuration variable is used as the new variable. The dual for quantum mechanical systems 
has been discussed i n We note however that there are important differences with the variables used for scalar 

field quantization in [121] . which do not use the available background metric in their definition. Our variables satisfy 
the canonical Poisson bracket 



U x (P^x))} = i\f{x)U x {P^(x)) . (5) 
A localized field may be defined by taking for example f(x) to be a Gaussian 

G(x,x fe ,/3) = e-^-^ 2 , (6) 

which is sharply peaked (1 <C 1 ) at a point Xk- We assume this in the following and write (ficixk) = 4>k- 

Since the representation of these variables that we will use is based on that for polymer quantum mechanics [l8l . [l|| , 
we briefly review this quantization before generalizing it to field theory. The Hilbert space is the space of almost 
periodic functions, where a wave function is written as the linear combination 



N N 
k=l k=l 



(7) 
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Here, the set of points {x k } is a selection (graph) from the real line. The inner product is 

_ r T 

T-,00 2T 



x\x')=\bn±-l dp e-**^' = 6 X , X ,, (8) 



in which plane waves are normalizable (the right hand side is the generalization of the Kronecker delta to continuous 
indices). On this Hilbcrt space the configuration and translation operators act as 

_ i^L e ipx* = Xk e lvx \ ^e ipXk = e lp{xk+x) . (9) 
op 

The momentum operator is not defined on this Hilbcrt space. Consequently, only the finite translation can be realized 
rather than the infinitesimal translation. 

The generalization of this quantization to field theory is straightforward. The Poisson bracket of the basic scalar 
field variables is realized as an operator relation on a Hilbert space with basis states 

\ai,a 2 , ■ ■ ■ ,a n ) (10) 

where the real numbers ai represent scalar field values at the spatial lattice points Xi. The inner product is 

(a'i,<4> ' ' ' ,a' n \ai,a 2 , ■ ■ ■ ,a n ) = S a ^ ai ■ ■ -5 a ^ an (11) 

when two states are associated with the same set of lattice points; if not the inner product is zero. For our purpose it 
suffices to consider a fixed spatial lattice of the type used in numerical computation. This inner product is background 
(metric) independent in the same way as for example that for Ising model spins; the difference is that for the latter 
there is a finite dimensional spin space at each lattice point, rather than a real number representing a field value at 
each point. In contrast the usual Fock space quantization uses the metric dependent Klein-Gordon inner product. 
The configuration and translation operators are defined by the following action: 

4> f \ai,a 2l ■ ■ ■ On) = aif(xi)\at,a2, ■■■a n ), (12) 

i 

U\(P^(xk))\ai,a 2 , ■••o n ) = |oi,02, ■•■ ,a k + A, • • -a n ). (13) 

It is readily verified that the commutator of these operators is a faithful realization of the corresponding Poisson 
bracket. The parameter A represents the discreteness scale in field configuration space. In this representation the 
momentum operator does not exist. There is however an alternative A dependent definition of 'momentum' given by 

P ^) = ^ [Ux(x k ) - Ul(x k )] , (14) 
which can be used to define the Hamiltonian operator. 



B. Semi-classical states 



Our goal is to obtain an effective Hamiltonian using this quantization procedure. This will be done by using 
semi-classical Gaussian states peaked at a classical phase space configuration. Such states have been defined for 
Friedmann- Robertson- Walker cosmology in the same representation |j20j , and we use those results here with a suitable 
generalization to field theory. 

For definiteness, we consider background spacetime here to be Minkowski spacetime with ^Jq = 1. This is also the 
focus of the present paper. In this spacetime consider a space lattice point x k and the basis states \a k = m\) k at this 
point, where m is an integer. A semi-classical state at x k is defined by 

00 

l^o>o>^ = ^ E e-^ 2 / 2 ° 2 e mX *°e- imXP °\m\) Xk . (15) 

m=— 00 

This is a Gaussian state of width a where the (real) parameters Po and </>o represent field values corresponding to a 
classical configuration at the point Xk- The normalization constant C > is given by the convergent sum 

oo 

C 2 = e- x2m2 ' a \ 2 ^ m . (16) 
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Calculation with this state gives the following expectation value [20J 

(Ux) = e lXPo{xk) e - x2/4a2 K(X,a,4> ) 

where 

K(\,a,(f>o) -- 



1 + 2 E™#0 COS 


A 


V 1 + 20 o <r> j 




e -7T 2 m 3 a- 2 /A 2 


1 + 2 Em^O C0S 


A 


e w 


m 2 er 2 /A 2 



Equation (fT7|) together with the definition (fl4|) gives the expectation value 

/ r>\r ^^ sin[P (x fc )A] _ A 2 /4(T 2 



This formula has the limits 



A 



lim {PUx k )) = sin[P (x fc )A]/A, 
lim (P A (a; fc )) = P (x fc ). 

A — »0 



(17) 
(18) 

(19) 

(20) 
(21) 



The first limit i.e. ignoring the width corrections of the quantum states, shows that the semi-classical state on the 
field lattice is peaked at the corresponding phase space value. The second limits shows that in the field continuum 
limit of the momentum expectation value has the appropriate classical value. This is crucial for the quantization 
scheme to be a viable quantization method given only finite field translation operators exist in the representation we 
are using. 



C. Effective Hamiltonian 



We now turn to derive an effective Hamiltonian for the scalar field using the semi-classical states (fT5|) . Let 
H{P < p(x),4>(x)) be the Hamiltonian density and we define the effective Hamiltonian density Ti cS by the expecta- 
tion value 

n eS (P (x k ),Mxk);a,X) = ~ (P o ,0o|(n A ) 2 +V^\P ,4>o)Z X ■ ( 22 ) 
The kinetic part of the operator gives 

1 



Po 



'o,MP A ) 2 |P o >o)SL A = y ( P o^o\ [2 - U 2X (P^x k )) - Ul x (P^(x k ))) |Po,<M^ A 



' sin 2 (PqA) +0 (—) , (23) 



A 2 



where the last term represents the width corrections of the semi-classical states. The scalar gradient operator can be 
defined by an expression such as 



dx 



1 - 

-<Px k +e 



(24) 



where cj> Xk is defined with a Gaussian peaked at x k as described above, and e is the spacing of a uniform lattice. (Note 
that there are two lattices in this formulation, one in field space with spacing A and the other in physical 3— space 
with spacing e.) Thus for wavelengths large compared to e the space continuum approximation applies, and we can 
write 

H cff (x) = ^sin 2 (PA) + i V 2 0, (25) 

where we have dropped the subscripts on Po and 4>o- I n this expression and what follows, we have neglected corrections 
due to the finite width of the semi-classical state. That is, we will work in the a — > oo limit for the rest of this paper. 
Finally, the effective Hamiltonian 

H cS = J d 3 x H cff (P,0) (26) 
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lives the modified wave equation 



d 2 4> 

dt 2 



1 - 4Af 



50 
~dt 



1/2 



(27) 



where we have set A = A* as a fixed "polymerization scale" . The nonlinear term represents the continuum remnant 
of polymer quantization. This equation is not Lorentz invariant (due to the (f> 2 term), but is invariant under Galilean 
and parity (x — > —x, t — > —t) transformations, as well as field reflection, (f> — > —<f>. We note also that the effective 
Hamiltonian is the integral of local density so local casual evolution via Hamilton's equations is assured. Indeed the 
corrections to the Hamiltonian in an expansion in A* gives terms that are higher derivative in momenta, which are 
like higher derivative terms from the lagrangian perspective. 

A generalization of this equation to arbitrary space metric q a b with lapse N = 1 and shift N a — derived using 
the same considerations appears in Appendix [A] 



D. Homogeneous scalar field on a cosmological background 

The formalism of polymer quantization presented here can be applied to different background spacetimes. In 
particular, it is useful to see what happens if we consider a reduction of the theory to the homogeneous spacetimc 
where it is a quantum mechanical system. Consider the Friedmann-Robertson- Walker (FRW) background 

ds 2 = -dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ), (28) 

so that y/q = a 3 where a is the scale factor. The basic observables (j4|) therefore reduce to 

c(> f = Voa 3 <j>, U x = e lXp l a \ (29) 

where Vq is a fiducial volume obtained from the spatial integration, and we have set / = 1 (which is natural given that 
any other non-zero constant can be absorbed in Vo- The algebra obtained from the fundamental bracket ([5]) gives 

{<p,U x } = -^U x . (30) 

These relations have an interesting feature when compared to the so called "improved dynamics" plj scenario 
in loop quantum cosmology (LQC) |22| : The appearance of the a 3 factor in U\ and the reduced Poison bracket is 
natural from the start, given the definition of the basic variables that we are using. This feature leads to the so-called 
dynamical lattice refinement in cosmology that is necessary to reproduce correct classical limit. In the LQC context 
on the other hand, there is no natural way to see how these factors can arise from the full theory upon reduction of 
a holonomy, which has no metric (y/q) dependent factors in it. 

Using these variables one can study the effective dynamics of a homogeneous polymer scalar field coupled to an FRW 
background, with rather surprising results: there is singularity avoidance and a built-in inflationary phase, without 
any need to quantize gravity [23j ] . In this scenario, the non-singular universe is past inflating. This is in contrast to 
LQC where the singularity avoidance is realized through a quantum bounce that occurs when the density of matter 
reaches the Planck density [2l|, |2J, |25|. Consideration of inhomogeneous fluctuations in this inflating background 
can provide a probe to study the polymer effects through the power spectrum of cosmological fluctuations. This is 
presently under investigation. 



III. PERTURB ATIVE SOLUTIONS 



A. General solution 



Our goal in this section is to develop perturbative solutions to the equation (|27[) . As a first step, we rewrite the 
equation in terms of dimensionless coordinates and fields. Noting that the parameter A* has dimensions of length 
squared, we define a mass scale M* by 



(31) 
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Using the mass scale M*, we can define a set of dimensionless quantities as follows 

$ = 4>/M k , T = M*t, X i = M k x i . (32) 

In terms of these variables, the equation (f27f becomes 

$- (1 -4$ 2 ) 1/2 D 2 $ = 0, (33) 

where an overdot indicates d/dT and D; = d/dX l . Notice that neither A* nor M* appears in this equation. The 
nonlinear wave equation (|33|) admits the trivial solution <f> = 0. We construct a perturbative series about this exact 
solution as follows: 

oo 

$ = e$ 1 +e 2 $ 2 + --- = ^e%, (34) 
i=l 

where e <g; 1 is a small dimensionless parameter. If we put this ansatz into the wave equation, expand in powers of e, 
and set the coefficient of each power equal to zero, we obtain: 

□ $i = 0, (35a) 

□ $ 2 = 0, (35b) 
□$ 3 = 2$ 2 D 2 $i, (35c) 

□ $ 4 = 2$i($iD 2 $ 2 + 2<j> 2 D 2 $i), (35d) 

□ $ 5 = 2($t + $ 2 + 2$!$ 3 )D 2 $i + 2$i($iD 2 $ 3 + 2$ 2 D 2 $ 2 ), (35e) 



where □ = -<9f, + D 2 . 

Before moving on, a few comments about these perturbative solutions are in order. First, note that the perturbative 
equations of motion (|35|) imply that we may absorb the $2 contribution into $1 by making the substitution $i+e$2 — ► 
$1. That is, we can effectively set $2 = 0. This then implies that we can absorb the $4 into $3 via $3 = $ 3 + ei>4, 
and so on. Thus at the end, we are left with $ 2 i = and a perturbative series with only odd powers of e. 

Secondly, it is perhaps useful to elaborate on the nature of our perturbative expansion. In field theory, perturbative 
series solutions often make use of dimensionless factors appearing in the action or equations of motion as expansion 
parameters. The prototypical example of this is a scalar field with a V((f>) — /i</> 4 self-interaction. In such a model, 
it is natural to develop a perturbative solution in terms of the dimensionless parameter /i; i.e., <p = ^2 n ^ n 4>n- Our 
situation is somewhat different since the wave equation (|33[) has no small dimensionless parameters appearing in 
it. Hence, we instead develop a series solution in terms of the amplitude of <I>. Such expansions are familiar from 
elementary physics; for example, consider the equation of motion of a simple pendulum of length £ in a uniform 
gravitational field g: 

(9= -sin (9, (36) 

where the dimensionless time is T = (jg/l) 1 ' i t = u;t and 9 measures the angular position. Like (f3"3")h this equation has 
no naturally small parameters appearing in it. To solve for 6, we employ an expansion 

00 

0(T) = £V0 n (T), (37) 

71=1 

which leads to equations of motion 

6q + 01=O, 03+02 = 0, 03 + 03 = g0?, (38) 

and so on. These are very similar in structure to (|35[) . and we see that 9\ satisfies the familiar equation of motion for 
a simple harmonic oscillator. The parameter e has a natural interpretation when we write down the full solution 

6 = ecas(ujt + 6) +■■■ (39) 



That is, e is just the dimensionless amplitude of the leading order term, which implies that the perturbative series 
is only valid for small angle deviations from the vertical. In our perturbative solution of (|33p we have precisely the 



7 



same interpretation: e is the amplitude of the leading order solution for $ and our perturbative series is only valid 
for small amplitude fluctuations. 

Finally, we mention that to re-introduce dimensions correctly one must multiplying (|34|1 by M+: 



£ 

i=l 



2i+l- 



(40) 



Note that <f> K<3 satisfies the dimensionful Klein-Gordon equation (— df + V 2 )(j) KG = 0. 

B. Planar symmetry 

We now seek solutions of the perturbative equations of motion (|35p with planar symmetry. That is, if we write 
the dimcnsionless spatial coordinates as X % = (X,Y, Z), then all the fields depend on T and X only. Defining the 
retarded and advanced times as 



U = T — X, V = T + X, 



(41) 



respectively, we assume that the leading order contribution $i is a plane-fronted wave packet propagating in the 
positive X direction: 



$i = F{T-X) = F{U) . 



(42) 



Here, F is an arbitrary function. As mentioned above, we can take $2 = without loss of generality. For $3, we find 
that (I35cl) reduces to 



with general solution 



$3 = 9l(V) + g R (U) - -V[F'(U)} 3 , 
6 



(43) 



(44) 



where gi, and c/r are arbitrary functions. To fix these functions, we need to choose initial conditions. A physically 
interesting situation involves preparing the total field $ such that the polymer corrections are zero on some initial 
hypersurface T,q; i.e., $ corresponds to a wave packet that initially evolves according to the Klein-Gordon equation. 
Obviously, the particular choice of the initial data hypersurface will affect the form of the solution. For simplicity, we 
choose Eq to be the null hypersurface V = and we require 



$3 =0 
^ 3 1 v=o 



This means that <7l(V) — and gn.(U) = 0. Restoring dimensions, we obtain the following solution for 



v M 2 fdf 



6/ M£ \du 



(45) 



(46) 



where 



M = eM*, u = t-x, v = t + x, f = S(u) = F(M*u). 



(47) 



Here, M is a mass scale representing the amplitude of the leading order (classical) solution. We now consider a couple 
of explicit solutions which are of physical interest. If the leading order term in the series (|4"6")) represents a Gaussian 
wave packet / = e~@ u e lku where the width of the wave-packet is 1//3, we obtain 



Me 



-(3 2 u 2 ^iku 



6 \" k 

Conversely, if the leading term represents a plane wave / = e lhu — e lfe (* -x ) then we get 

1 M 2 k 2 



= Me lHt ~ x) 



1 + r 



6 Mf 



-ik(t + x)e 2lki - t ~ x) 



(48) 



(49) 
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In either case, we see that the polymer corrections are small if 

kv=k(t+ X )«^y ^ 2 . (50) 

Naively, we recover Klein-Gordon solutions in the limit of small amplitudes M <C A/* and/or long wavelengths k <C M*. 
However, notice that when (t + x) becomes large, the leading order polymer correction grows to become comparable 
to the Klein-Gordon contribution. In other words, the perturbative solutions will eventually break down. This 
implies that if we wait long enough, then the nonlinear quantum effects will affect the wave propagation significantly, 
irrespective of the values of k or M. 

Previously, we considered the initial data hypersurface to be null. We now consider the situation where initial data 
are specified at the T = hypersurface. As earlier, we require "purely Klein-Gordon" data on the initial hypersurface. 
That is, <I>3 should satisfy the following conditions 

$ 3| T=0 = 0, $ 3 | T=0 -0. (51) 

Here it is slightly more involved to find explicit expressions for <7l and g^ . In the case where the leading order term 
represents a plane wave, we get 



<t> = Me lk{t -' x) 



1 + e-««8in3fct 



3 A/* V 3fci 

It follows from the expression ([52")) of perturbative solutions that polymer corrections are small only if 



(52) 



(»»> 



M J \ k 

This is effectively the same criterion as found earlier for the initial condition specified at null hypersurface (|45|) . 
In particular, we see that after a sufficiently long time the polymer corrections will become non-negligible and the 
perturbative solutions will break down. 



C. Spherical symmetry 



It is also interesting to solve the perturbation equations (|35[) under the assumption of spherical symmetry; i.e., the 
fields are functions of t and a radial coordinate r. If we assume that the leading order contribution represents an 
outgoing spherical wave packet, we obtain the following solution: 



Me -p 2 u 2 e iku 
kr 



2i k 2 M 2 
~3~ Mi 



I 

kv 



1 

kv 



2i(3 2 u 



e -2/3^iT e 2iku 



(54) 



t — r and v — t + r. We have imposed the boundary conditions such that <j> reduces to the Klein-Gordon 



Here, u - 
solution . 
the solution reduces to 



Me /kr on an initial null hypersurface v = vq. It is interesting to note that for late times (i.e, v — > oo) 



Me -p u e iku 
kr 



1 



2i kM 
~3v Mi 



I 



2if3 2 u 



e ' e 



(55) 



which is actually a solution of the ordinary Klein-Gordon equation (— d 2 + V 2 )0 = 0. In contrast to case of planar 
waves, the polymer effects on spherical waves gradually decreases as the wave packet travels away from the initial 
surface. The reason behind this is that the effective amplitude of the Klein-Gordon part of the solution decays as 1/v 
when it is far away from r = 0. Since the polymer corrections are explicitly amplitude dependent, it follows that we 
recover Klein-Gordon behaviour as v — > oo. It may be emphasized that this behaviour however does not mean that 
the asymptotic form of <p exactly matches the Klein-Gordon prediction. In particular, the polymer effects generate 
a whole series of additional wave packets with frequency 3fc, 5fe, etc in addition to the original wave packet with 
frequency k, The relative amplitude of the "polymer harmonic" with k n = (2n + l)k will be roughly (kM 2 /vqM^) 12 , 
where Vq can be taken as the initial radius of the wave packet. From this, it follows that our perturbative solutions 
will be valid at infinity if 



Ml « { ( M 



(56) 



That is, low wavenumber and low amplitude waves are well described by perturbative techniques. 
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IV. PHENOMENOLOGICAL IMPLICATIONS 



In this section, we explore some phenomenological implications of the perturbative solutions of the semi-classical 
Klein-Gordon equation. We begin by analyzing the propagation of a planar wave-packet and derive its group velocity 
and an effective "dispersion relation" (which must be interpreted with care). We then examine the conditions under 
which our perturbative solutions break-down, and hence derive a "coherence length" over which planar waves resemble 
solutions of the ordinary Klein-Gordon equation. Finally, we consider the spherical wave-packet solutions in the context 
of polymer modifications to the spectrum of gamma ray bursts. 



A. Group velocity of polymer wave-packets and effective dispersion relations 

As mentioned earlier, there has been much discussion in the literature concerning "modified dispersion relations" 
of the form |T]), especially in the context of doubly special relativity [26|, [27|]. Such formulae are consistent with 
linear wave equations with higher derivatives; for example, an equation of the form [— d 2 + V 2 — ^piV 4 ]</> = leads 
to the dispersion relation lu 2 = fc 2 (l + ^i^ 2 )- However, for nonlinear wave equations such as the polymer-corrected 
Klein-Gordon equation (|27p . the derivation and interpretation of dispersion relations is more complicated. In this 
section, we suggest a physically-motivated derivation of an effective dispersion relation based on the propagation of 
spatially localized wave-packets. 

Let us begin by considering the polymer corrected planar wave solution (|46p . We assume that the leading order 
Klein-Gordon term in the series represents a spatially localized waveform with characteristic frequency k: 

f(u) = h{ku)e lku . (57) 

The real and dimensionless function h defines the envelope of the wave-packet. For simplicity, we restrict ourselves 
to symmetric wave-packets with h(£) — h(—£) and define M such that h(0) = 1. Also, we assume that h vanishes 
strongly at infinity such that integrals of the form J £™/i 2 (<!;)<i£ are finite. Sufficient conditions for this hold are that 
h has compact supports or decays exponentially for large arguments. Note that our choice of f(u) means that the 
Klein-Gordon contribution to the perturbation series </> K g = 

Mh(kuy ku is self-similar; i.e., KG is invariant under a 
rescaling of the coordinate u — > £u if we also rescale the wavenumber k — > k/L 

In this subsection, we are interested in the leading order polymer effects, so we truncate the series (05]) as follows: 



<j>^Mh{ku)e^- -—^—[h(ku)e iku ]j . (58) 
We define the position of the packet at a given time t by the average value of x: 

_ JZc dx x _^± _, $duu<j>*<j> 

where we have made use of u — t — x. In general, x avg will depend on time and the wavenumber k of the packet. It 
is easy to confirm that in the Klein-Gordon limit M* -^oowe get x avg = t; i.e. the pulse propagates at the speed of 
light. 

To calculate x avg in general, it useful to define 

Q(0 = h(0e*{^[h(0e-*]\ 3 + c.c. (60) 
From the fact that h is even, it follows that Q is and odd function of £. We also define the dimensionless moments 



/'r, 



/oo /*oo 
d£ r* 2 (0. "n = / d£ CQ(0 , (61) 



with n = 0, 1, 2 . . . From the parity of h and Q, these satisfy M2n+i — and v 2n — 0. Using these definitions and the 
fact that v = 2t — u, we find that 

#* = M 2 \h(ku)\ 2 k3{2t - u) (^) 4 Q(M + O (i) , (62) 
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M 2 t ( M\ 4 1 / M \ 4 „ / 1 



which leads to 

This gives the effective position of the pulse as a function of time 

M 2 k 2 vi\ m ( 1 



^=T + w^ +0 UJ- (64) 

From this expression, we see that the polymer wave-packet travels with a group velocity different from unity 

dx mg M 2 k 2 / 1 \ _ v x 



dt M$ \M* J ' 3no 

Here, Th is a dimensionless form factor that depends on the wave-packet profile h only. For example, if h is a Gaussian 
we obtain 

HO = <r ei "' * *-£(2£2)«-». (66) 

We see from f|65[) that the polymer correction to the group velocity depends not only of the wavenumber k, but also 
on the Klein-Gordon amplitude M and the shape of the pulse Th- This formula for v g can be physically interpreted 
by noting that the characteristic energy density associated with the Klein-Gordon part of the wave-packet is 

Pkg = \ [(^kg) 2 + (V0 KG ) 2 ] ~ M 2 k 2 . (67) 

Then we see that 



(68) 



That is, the polymer correction to the group velocity scales like the density of the wave-packet normalized by the 
"polymer density" Mf . This is somewhat akin to the situation in loop quantum cosmology, where polymer corrections 
to the Friedmann equation scale like the density normalized by the gravitational polymer scale M^. Previously in 
the literature, the astrophysical and experimental consequences of modified group velocities have been derived by 
considering corrections that depend only on k Q- This assumption was motivated by ad hoc dispersion relations 
of the form JT]). Here, we see that the corrections can depend on the amplitude and shape of the wave-packet as 
well. This crucial difference is due to the fact that we are dealing with a nonlinear semi-classical wave equation. It is 
worthwhile stressing that the above group velocity has been derived, not postulated as in previous work. 

Given the above expression for the group velocity, one may be tempted to derive an effective dispersion relation by 
setting v g = dio/dk and integrating with respect to k. This can be done, and we find 



,2 



k 



2M 2 k 2 / 1 



(69) 



On first glance, this seems to be of the same form as ([T]) with 7 = 2. But we again note that corrections depend 
on the amplitude M and shape Th of the original wave-packet. We comment that amplitude dependent dispersion 
relations are a common consequence of nonlinear wave equations; for example, they arise in the theory of deep water 
surface waves [28| . 

We conclude by noting that significant caution is warranted if we try to interpret (|69[) as the dispersion relation for 
some linear wave equation with higher derivative terms. For example, the fact that d 2 uj/dk 2 ^ would naively imply 
that our wave-packet is spreading as time progresses. But an explicit calculation of the effective width of the pulse 
gives 



AX = -^ = ^ + 6M^UJ V 3 -—) +0 {W?)- (70) 

Note this is time independent to leading order in M+. So although the dispersion relation insinuates that the packet 
is spreading, a direct calculation shows that this is not the case to leading order. Fundamentally, simple dispersion 
relations of the form ([T]) are not capable of capturing all of the interesting dynamics associated with our nonlinear 
wave equation (|2"7| . 
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B. Effective coherence length for planar waves 



In this subsection, we discuss the break-down of the perturbative solutions to the polymerized wave equation and the 
associated observational effects. As we have seen in Sections IIII Bl and IIII CI the polymer quantum effects accumulate 
or dissipate with time for planar and spherical waves, respectively. So, to maximize our chances of seeing deviations 
from conventional theory we should look at highly collimated beams of radiation. To quantify how big polymer effects 
are in such beams, we rewrite our polymer-corrected plane wave solution (|49[) as 



2poly 



2ik(t-x) 



where 



= Me'*'*" 1 ' 



l/A^y/M^ 2 

dpoly - ipoly - k I M J ^ k 



(71) 



(72) 



As we have mentioned before, the magnitude of the polymer corrections inside the square brackets becomes larger as 
time goes on, and our perturbative solutions are valid for 



t <§C ^poly 



(73) 



That is, when this condition is satisfied the approximation cj> <j} KG holds. We note that we can obtain a very similar 
expression to (|71|) if (f> KG is a wave packet whose width is much larger than 1/k. In that case, the validity of the 
classical solution demands that the distance d traveled by the wave packet satisfy 



d < d 



poly • 



(74) 



It is useful to think of t po iy and d po iy as a "coherence time" and "coherence length" , respectively, because they 
represent the time or distance over which the polymer corrections to a classical wave packet remain small. 

More practical expressions for i po iy and d po i y come from noting that the luminosity (particles per unit area per unit 
time) associated with the Klein-Gordon part of the beam is 



Pkg 

k 



1 

2k 



[(cWkg) 2 + (V0kg) 2 ] ~M 2 k 



Restoring conventional CGS units, this yields 

4 / x -2 

/ \ / \ 

^poly " ^pl 



Mpi 



E K 



E, 



pi 



L p i — 



7 x 10 108 cm- 2 sec _1 



(75) 



(76) 



where E KG = Hk/c is the energy per particle in the beam and L„\ is the Planck luminosity. From this, we see that the 
coherence length is minimized for high-energy and large luminosity beams. 

One of the most luminous artificial beams ever constructed is in the Large Hadron Collider (LHC) 29]. In that 
device, the proton luminosity is L KG ~ 10 34 cm -2 sec -1 at an energy of E KG ~ 7TeV. This corresponds to 



Ipoly 



pi 



10 45 Gpc . 



(77) 



The coherence length is longer if we consider the most energetic cosmic rays: The Pierre Auger Observatory observed 
561 particles of energy ~ 10 10 GeV with an integrated exposure of ~ 10 5 km 2 yr [3C| . which corresponds to 



Ipoly 



10 84 Gpc 



(78) 



Unless M* <C M p \, both of these distances are many orders of magnitude larger than the size of the observable 
universe (which is ~ 3 Gpc). Hence, if the polymer scale is of the same order as the Planck scale, the Klein-Gordon 
equation gives an excellent description of the physics in both scenarios, and polymer effects will be very hard to 
observe. However, we should point out that there is no a priori reason to take ~ M v \ in this model. Indeed if 
we take = 5 TeV, we find that t po \ y — 1 sec in the LHC. That is, polymer effects would completely dominate the 
beam after only one second. Since we are unaware of any exotic beam behaviour in the LHC, this can serve as an 
effective lower bound: M* > 5 TeV. 
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C. Spectral modifications to gamma ray bursts 



We conclude this section by considering the spherical wave packets of Section IIII CI We can use these to model 
the radiation emitted from a gamma-ray burst (GRB) if we assume that the energy emitted by the progenitor is 
isotropic and we can approximate the electromagnetic wave dynamics using a scalar field. Both assumptions should 
be sufficient for the order of magnitude estimates we seek here. 

We imagine the following situation: The gamma-ray progenitor is located a distance of d away from the earth. At 
t = to a spherical pulse of radiation of radius ro = to is emitted from this progenitor. The initial pulse profile is 
assumed to be Gaussian: 



Me -f3 2 (t-r) 2 e ik(t-r) 

kr 



(t « t ,r w ro) 



Now, the signal S(t) seen by an earth based detector will just be the polymer solution (j54")) evaluated at r = d: 



S(t) 



\r=d 



Me- f) T e lkT 
kd 



i kM 2 



3 M*ro 



1 + 



2i/3 2 r 



-20 2 r 2 2ikr 



with t — t — d and assuming d S> ro . The Fourier transform of this signal is 



S(u) 



1 



dre-^S{r) = fte-C^W +iS a f-^-V 

\3k/ 



18(3 2 



-(w-3fe) 2 /12/3 2 



where 



Si = 



M 



V2f3kd' 



S,= 



V6 M 3 
1& pdr Mf 



(79) 



(80) 



(81) 



(82) 



The first term in (fBTj) represents the spectral line centered about u> = k that our detector would have seen without 
polymer quantization. The second term is the result of the leading order polymer correction: An additional spectral 
line centered about to = 3/c. The relative amplitude of the two lines is 



Si 



E tot k0 



r Q M$ ' 

where we have noted that the total energy in the pulse initially is 



(83) 



EL, 



Pkg d 3 x ~ M 2 /f3, 



(84) 



assuming a long duration k 3> (3 and a large initial radius ro/3 3> 1. Now, the initial pulse must have a size bigger 



than the Schwarzschild radius associated with E tot ; i.e., tq > E tot /M 2 x , so we obtain an upper bound on S3 /Si: 



S3 
Si 



< 



E K 



pi 



At 



M* 
M P i 



10 



-66 



E K 



At 



10 5 eV J V 1 sec 



(85) 



where At is the duration of the pulse observed on earth and E KG is the energy per particle (photon) in the beam; for 
a GRB, one typically has At ~ 1 sec and E KG ~ 10 5 eV [31]. As in the case for the LHC, we know of no evidence for 
exotic polymer effect in observed GRBs, which suggests that S3 /Si -C 1. This yields another bound: M* > ITeV. 
This is virtually the same numeric result we obtained from the LHC. 



V. DISCUSSION 



We have developed an approach for studying polymer quantum field theory on a curved background spacetime. The 
quantization procedure is motivated by ideas used in loop quantum gravity. However, the representation used here is 
significantly different to the one used earlier [T^]. The polymer quantization differs from the usual approach in two 
important aspects: the inner product is background independent, and there is a built-in mass (or length) scale M* due 
to the choice of basic observables. Using semi-classical methods we have derived an effective nonlinear wave equation 
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governing the field dynamics. The form of this semi-classical wave equation necessarily depends on the underlying 
quantum state, but we explicitly find that it reduces to the conventional Klein-Gordon equation at low energies. 

In §1111 we developed perturbative solutions to the nonlinear wave equation. The natural expansion parameter 
involves the amplitude of the lowest order contribution. We then specialized to waves with planar and spherical 
symmetry in mil Bl and mil CI respectively. For the plane-symmetric case, we found that the polymer corrections 
to the classical solution actually grow in time. This suggests that if one waits long enough, nonlinear effects will 
completely dominate the Klein-Gordon contribution to the wave. 

The converse is true for spherical waves, which tend to decrease in amplitude as they travel outward from the center 
of symmetry. This implies that the polymer effects remain bounded as time progresses. This qualitative difference 
between planar and spherical wave highlights one of the most important features of the model: corrections to classical 
dynamics depend on more than just the frequency of the wave. Furthermore, if one begins with a spatially localized 
spherical wave-packet with characteristic wavenumber k, one obtains a whole series of wave-packets with wavenumber 
k, 3k, 5k . . . at infinity. 

In £HV( we derived some phenomenological implications of the effective wave equation. We found that the group 
velocity of planar wave-packets depends on their frequency, amplitude and shape. This is in sharp contrast to many 
models in the literature, which assume that the group velocity only recieves frequency-dependent modifications. From 
the these group velocities, we derive effective dispersion relations for the model. These bear some resemblance to the 
dispersion relations of weakly nonlinear waters waves in that they depend explicitly on amplitude. Our method of 
extracting a dispersion relation appears to be new, and has potential applications to other nonlinear wave equations. 

The fact that our perturbative solution breaks down after a finite amount of time has interesting practical im- 
plications, as discussed in SjIVBj This implies that there is a "polymer coherence length" c? po i y associated with 
plane-symmetric wave packets. After a wave packet has traveled a distance d po i y , polymer effects will completely 
dominate the classical behaviour. We express d po i y as a function of the frequency and luminosity of the beam. We 
calculate G? po i y for the proton-antiproton beam in the Large Hadron Collider, and use the fact that no exotic behaviour 
has been reported in that device to deduce a lower limit of M* > 1 TeV. We also calculated d po iy for a highly colli- 
mated beam of cosmic rays, and found that it is many of orders of magnitude larger than the size of the observable 
universe. Finally, we considered spherical gamma ray bursts in fllV CI and demonstrated that if the characteristic 
frequency of the burst is u> = k, an earth-based observer will see additional spectral lines centered at lu = (2n + l)k. 
From the amplitude ratio of the classical to polymer spectral features, we deduced that > 1 TeV consistent with 
the LHC result mentioned above. 

It is interesting to note that all of the modifications to wave propagation mentioned above are derived without 
including any quantum gravity effects, unlike other approaches such as doubly special relativity, loop quantum gravity, 
or non-commutative spacetime models. Another important feature worth emphasizing is that fundamental discreteness 
in this model leads to nonlinearities at short distances. This is contrary to the common assumption that the effects 
of fundamental discreteness manifest themselves via a linear wave equation with higher derivative terms. 

There are a number of possible developments based on the ideas we have discussed, such as applications to Hawking 
radiation and to the spectrum of fluctuations in inflationary cosmology, which have as their basis a wave equation on a 
curved background. These are presently being studied. Another application is to homogeneous scalar field propagation 
on an FRW background. This work provides some novel results, including cosmological singularity avoidance and 
a mechanism for inflation [23j . The former result is a concrete realization of the intuition that in the Hamiltonian 
constraint, the compactification of the matter energy density that results from the semi-classical effective Hamiltonian 
leads to a bound on the gravitational terms in this constraint; the converse occurs in the LQC case, where it is the 
gravitational kinetic terms that are compactified by polymer geometry. We also show in [23| that the energy scale of 
polymer inflation is fixed by M+, which leads to an entirely different way of constraining the fundamental discreteness 
scale in our model. 



Acknowledgments 

This work was supported in part by the Natural Science and Engineering Research Council of Canada, and the 
Atlantic Association for Research in the Mathematical Sciences. 



APPENDIX A: THE SEMI-CLASSICAL WAVE EQUATION ON A CURVED BACKGROUND 

The method followed to derive an effective wave equation on Minkowski spacetime is applicable also to an arbitrary 
background. Here we restrict attention to metrics with lapse N = 1 and shift N a = 0, and denote the spatial metric 
by q a b and its determinant by q. 
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Let Ti(P c f > (x), 4>{x)) be the Hamiltonian density and define the effective Hamiltonian density H cS by the expectation 
value 



\\2 



(Al) 



The semi-classical states given above have to be modified to take into into the density weights, and the fact that that 
the peaking value for the scalar field is actually for 0/ rather than <fr. With these changes, to leading order in the 
width cr, the kinetic part of the operator gives 



(Po,M^ A ) 2 |Po,0o)SL A = £ (Po,M (2-U 2X (P4x k ))-Ul x (P^x k ))) |Po,0o)^ 



q . 2 W 



oi 1 - 



(A2) 



where the last term represents the width corrections of the semi-classical states. The scalar gradient operator is 
diagonal in the basis, and as for the flat space case, its discrete aspect can be ignored for wavelengths large compared 
to the spatial lattice spacing e. Thus we can write 



n/ctt/ \ yl ■ 2 



PA 

71 



(A3) 



This Hamiltonian density leads to the modified wave equation 



d 2 (j) 



1-AX 2 





1/2 







4Ag \ dt 



arcsin 2 A 



' dt 



= 0. 



(A4) 
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